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Abstract 

In this paper we study the non-Gaussian features of the primordial fluctuations in loop quan- 
tum cosmology with the inverse volume corrections. The detailed analysis is performed in the 
single field slow-roll inflationary models. However, our results reflect the universal characteristics 
of bispectrum in loop quantum cosmology. The main corrections to the scalar bispectrum come 
from two aspects: one is the modifications to the standard Bunch-Davies vacuum, the other is the 
corrections to the background dependent variables, such as slow-roll parameters. Our calculations 
show that the loop quantum corrections make fy, of the inflationary models increase 0.1%. More- 
over, we find that two new shapes of non-Gaussian signal arise, which we name 71 and 7». The 
former gives a unique loop quantum feature which is less correlated with the local, equilateral and 


single types, while the latter is highly correlated with the local one. 
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I. INTRODUCTION 


As a non-perturbative and background-independent theory, Loop Quantum Gravity 
(LQG) has achieved great successes in past years: derivations of the quantized area and 
volume operators , calculations of black holes entropy d and Loop Quantum Cosmol- 
ogy(LQC) |9], etc. And the nonperturbative quantization procedure of LQG is also valid for 
a more general class of four-dimensional metric theories of gravity OU As an example 
of LQG, LQC gives a quantization scheme of LQG for a symmetry-reduced model in the 
homogeneous and isotropic Friedmann-Lemaitre-Robertson-Walker universe. The discrete 
spacetime geometry in LQC scenarios predicts a non-singular bouncing universe in some 
simplified models, which satisfies most of the astronomical and cosmological observational 
constraints. Although the quantum correction effects are being diluted with the expansion 
of our universe, it remains present in a weaker form, especially on/near the super-horizon 


scales. 


for inverse volume corrections and in |16 


Recently, the gauge invariant cosmological perturbation theory has been systematically 
constructed in 


17] for holonomy correc- 
tions. Some relevant applications have been considered in |18,|19}, |21),|22]. The inverse volume 
corrections and the holonomy corrections are two main quantum corrections in LQC. The 
inverse volume corrections come from the quantization of the inverse of the volume operator 
in LQG. The inverse volumes exist in the Hamiltonian constraint of gravity and the usual 
matter Hamiltonian, especially in the kinematic terms. Since the volume can be taken the 
value zero, there does not exist well defined inverses of the volume operator. Fortunately, 
with the Thiemann trick , we can construct well defined inverse volume operators, which 
bring the quantum corrections. While the holonomy corrections arise from the loop quanti- 
zation based on the holonomies instead of the direct connection. The holonomy corrections 
become important when the energy scale of our Universe approaches the Planck one. Both 
the modifications to the scalar and tensor primordial power spectra from the inverse volume 
corrections are carefully investigated by the authors of [18]. Their results show that the in- 
verse volume corrections could give rise to the enhancement of the power spectra on the large 
scales, i.e., a red-tilt one. However, some other mechanisms such as the non-commutative 


geometry or string theory lod bd bd. could also lead to similar features. Therefore to seek 


for signature of loop quantum cosmology, the study for the non-Gaussianity features in loop 


201709.00169v1 


chinaXiv 


quantum cosmology is necessary. 

Because the primordial non-Gaussianities are quite helpful to distinguish inflationary 
models, so far a lot of papers have been devoted to studying the non-Gaussianities in different 
inflation models, see the relevant references in the nice reviews bd . Inspired by the 
studies of "n in this paper we mainly consider the non-Gaussianities from the inverse 
volume corrections in LQC. The reason for consider inverse volume corrections only is as 
following. We denote ðiny as the correction term coming from the inverse volume operator 


and änt as the correction term from the holonomy corrections. We can estimate the inverse 


Har 2 
Os deg (Zo ’ (1) 
3 pp 


where the Planck density pp; is assumed as the quantum gravity scale. From the above 


volume correction as|21 


expression, we can see that the inverse volume corrections behave very differently from what 
is normally expected for quantum gravity. For low densities, the holonomy corrections is 
small, but the inverse volume one may still be large because they are magnified by the inverse 
of dno. For an example, the small holonomy corrections of size dno, < 1076 then requires 
the inverse volume corrections larger than Au, > 10~° even at scale p ~ 10 ?pp,. These 
novel features make the investigations on the inverse volume corrections more interesting 
than the former at sub-Planckian inflationary scales. So we only consider the inverse volume 
correction in this work. 

Explicitly, in the perturbation theory in LQC, the inverse volume operator can be cap- 
tured by a correction function such as à ~ 1 + agói, ~ 1 + Q0(Giny/a)”, where a is the scale 
factor of the FLRW universe and a;,, is introduced to describe the characteristic scale of the 
inverse volume correction, which is not the Planck one in general. When aj,,/a < 1, we can 
ignore the correction term. However, if a;,,/a < 1 during inflation, one cannot neglect in- 
verse volume corrections. In this case, the correction approximates aodiny(k) 2 ó(ko)(ko/ k)^, 
where k and ko are, respectively, the considered perturbation wave number and some char- 
acteristic number involved in the inverse volume correction. In addition, many works about 
LQC imply that c € [0,6] [19]. From the form of the inverse volume correction, we can see 
that a small sigma corresponds to a small the inverse volume correct, vice versa. Therefore, 
as an example, following lid. we take o = 2 in this paper. For other c values the behavior 


will be similar. Furthermore, in terms of spherical multiples the wave number could be 
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expressed as k zz 10~*hl, with h representing for the reduced Hubble parameter h = 0.7 and 
l for the spherical multiples. In the typical linear regime of Cosmic Microwave Background 
(CMB), the multiples / range in 2 « | < 1000 or more. Given all the mentioned observations, 
we expect some new features in the non-Gaussianities will arise. The purpose of this work 
is to investigate the characterized sizes and shapes of bispectra in LQC scenarios. 

Note that the bispectra for the single field slow-roll inflationary model have been calcu- 
lated in the papers For simplicity and comparison, we study the simplest single 
field slow-roll inflationary model in LQC. Our results show that the quantum corrections 
mainly come from the third order interaction Hamiltonian, the corrected vacuum state and 
the corrections to the slow-roll parameters. Our paper is organized as follows. In Sec. 
the canonical formulism and the slow-roll inflationary model in LQC scenarios are briefly 
reviewed. In Sec. we study the power spectrum in LQC and recover the previous results. 
The effect of the inverse volume corrections on the non-Gaussianity in LQC is investigated 
in Sec. The detailed analysis for the sizes, shapes and shape correlations is presented 
in Sec. [Vhnd Sec. [V] respectively. Throughout this paper we set 87yG = 1 and Einstein’s 


summing convention is always adopted. 


II. REVIEW OF LOOP QUANTUM COSMOLOGY 


The framework of LQG/LQC will be briefly presented in this section. Firstly, we discuss 
the canonical formalism in LQG, and then introduce the dynamics of slow-roll inflationary 


models in LQC scenarios. 


A. The canonical formalism in loop quantum gravity 


In the framework of LQG H d. the spatial metric as a canonical field is replaced by the 
densitized triad E?, defined as 


E; = | det(e)|ef, (2) 


where e? is the inverse of the cotriad e, related to the spatial metric by qu, = efef. The 
canonically conjugate variable to the densitized triad is the Ashtekar-Barbero connection 


Ai := TI} + 4 Ki, where K? is the extrinsic curvature and y z 0.274 is the Barbero-Immirzi 
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— a bd. The densitized triad E" and the Ashtekar-Barbero connection A? satisfy 


the following commutator relation 
(Ale), Ej(y)) = 95050 (a, y). (3) 


The spin connection I? is defined such that it leaves the triad covariantly constant and has 


the explicit form 
T 1 
D, = EM e; (aly Sue + 66s Diei] Sm). (4) 


In the new Ashtekar variables, the Einstein-Hilbert action can be expressed in the canonical 


form 
Bos d Zi | d'a KE? — Gen [N*] — 2e [IN] — Geral] (5) 
S 
where the Diffeomorphism constraint is 
1 , 
SCHER n PrN’ JEE — à A E? — ABE). (6) 
by 
And correspondingly the Hamiltonian constraint can be expressed as 
NI ~ / Pend E Do Ti + d, (I"T? Seat (7) 
grav 2 e i [det E] cd mn\*~c~d c d/|* 
The Gaussian constraint is 
C NH I = | d? zA' (0E? + eT) Ez + e KIER), (8) 
x 


which can be solved through standard procedure A Thus, solutions for the scalar mode 
perturbations are completely determined by the Hamiltonian constraint and the Diffeomor- 


phism constraint. 


B. Slow-roll inflationary models 


In this subsection, we shortly review the inflationary dynamics of the Friedmann- 
Lemaitre-Robertson-Walker Universe in LQC scenarios. The modified Friedmann equation, 


Raychaudhuri equation and Klein-Gordon equation are respectively 


oc: 
w= ta(S +m). H 
142 aP) lo 2/4 Vp 
w = H? (1+ =) Dm (1 a): d) 
ei + wg! (1 Sall + DEA = 0, (11) 
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/ 


— P 
where H = 3p 


v, = dv/dp, p = a? and a prime represents the derivative with respect to 


the conformal time. à and Y are the correction functions for the inverse volumes and they 


read à z 1 + Ain pelt VoOinv 
Following lid b] the slow-roll parameter e can be straightforwardly calculated as 


H' 
€ 1— ag à 
= eo(1 sr ^feÓiny.) D (12) 


where eo denotes for the usual slow-roll parameter. And the explicit form of coupling constant 
^j, reads 


zoe [EB A 
"en dE E (13) 


Typically we set ou = 0.06, vo = 0.17 and e; = 0.01 in this paper lid. do is determined 
by quantum correction and the analysis after Eq.(43) implies that it takes O(10 ?). As 
mentioned in Introduction we set o € [0,6] in this work. From above expression we can 
easily estimate that ye is of the order O(10 ?). Formally we can also express the another 
slow-roll parameter 7 as 


n 


9-1-£ 
Hy!’ 
= ell Br Vn On) H (14) 


where no denotes for the usual slow-roll parameter as eo. 

As in the usual situation of the single field inflation model, we can assume that the two 
slow-roll parameters, €o and rp, take roughly the same order as 10~?. Thus the typical values 
of coupling constant 7, and ye are of the order O(10~*). The terms proportional to the diny 
represent for the inverse volume corrections. Here, we should emphasize that the subscript 


“inv” is introduced to avoid confusion with perturbations, such as dy. 


Ill. POWER SPECTRUM 


In this section, we will firstly review the formalism of scalar perturbations in LQC; then, 
derive the second order Hamiltonian; and finally calculate the primordial power spectrum 


in the spatially flat gauge. 
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A. Formalism on the scalar modes 


Consider the scalar perturbations only, the general form of a perturbed metric around 


the isotropic FRW background is 
ds? = e(74 — (1 + 26)dr? + 28, Bdrdz* + a — 2) ban + 20,065 | dz*da^| , — (15) 


where the scalar factor a is a function of the conformal time 7, and (¢, v, E, B) are the four 


scalar metric perturbations. In the perturbation theory, the triad can be described by 
E? = Ee + ôE? , (16) 
where 
En = 76%, SE? = —2pwo? + P(A — 0°0,)E. (17) 


The perturbed triad is described by the spatial part of the perturbed metric ai and E. Here 
A is the laplace operator in the flat space. Similarly, the perturbed lapse function and shift 


vector can be described by the other two scalar metric perturbation o and B respectively, 
AN = Nọ, ôN" = 0° B. (18) 
The extrinsic curvature can be perturbed as 
Ki = Ki + 8K! = béi 68K. (19) 
For a general triad (16), the linearized spin connection becomes 
srt = lE! 20 
a pe pU rj. ( ) 


As described above, the symplectic structure also splits into two parts, one for the back- 


ground variables and the other for the perturbations, 
fox 1 i i 
{k, p} SS 3V (Ak Loi, SE (y)} = S(x, y) 020; H (21) 
where the background variables are defined by 
p= — | Edr, k= — | Ree. (22) 


Here Vo is some artificial finite volume. 
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In this paper, the matter part is represented by a scalar field y. Similarly, we split the 
field y and its conjugate momentum 7 into homogeneous part and inhomogeneous one as 


well 
p=G+6p, m—m-cóm. (23) 
Hence, the basic Poisson brackets are reduced into 


{6,7} = c Iëetel, 6n(y)) = (x — y). (24) 


For simplicity, we introduce the LQC formalism with finite cell Vo rather than the whole 
R? region in above description. But the unphysical feature of Vo can be remedied by lattice 
refinement model . Since our following calculation only involves Aan, we adopt the lattice 
refinement parametrization procedure Ss to eliminate the effect of artificial volume Vo. 


B. The second order Hamiltonian 


According to , the quantum corrected second order Hamiltonian constraint can con- 


veniently be written as 


HO) = AECH IN] + 260 [5N] + Hr LN] + 260 aM 


grav grav matter matter 


1 - 1 V d 
2 [ PrNaH, + 5 J PLENAS, + f PaN v5 - 050 + a 
E D 2 


+ 


| PordN DH + $09] (25) 
M 
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` b D ` | 
$0. = —4(1 + f)kV/posóK — (1+ 9) GE; + "a ; (26) 
= j c = Ce 2k C j 
HO, = VpOKLSK GOOF — VPK)? — E 
ER c k Jy k? SEW 07 di c d 
e SE E + BP (GEI — (1+ haan zm (O.5E%) (846 Eg) , (27) 
TOT zi SES 
$0 = (Ui f)— EXE -+ fz a ; (28) 
$9 = 0. (29) 
Jo ES 
SQ = p (a= velee VIZ) | (30) 
On? ôn ðE? 17? (SIES)? AEOS ECS EA 
gaz Ge Ree quoe qc 
HY = (1 + nope (1 "qr g2 DI 2p 2p 2 p3/2 (c + 93) 8p? Ap? ` 
(31) 
1 
99 = zC + gs) vp? AS PASY (32) 
1 SSES (ISES)? AKSES EA 
—3/2 HE 2 ci c = c Cd k 
Hy) = STT IO g) 5 Viool PS + Volpi + V(?) [LI am p 
(33) 


where the definitions of the counterterms can be found in or in the Appendix B of i4. 


And the perturbed second order diffeomorphism constraint can be expressed as 


gP WAN" = go AN" + 20. WV" 


grav 


d d*xóN^ [pos (863) — p(O,6K*) — kó*(O48E2) + (18,3) .— (34) 
X 


Based on this corrected Hamiltonian, we can get the homogeneous and inhomogeneous part 


of matter field as 


o ETA 
äise =p} [ie am foe] o) +p y (35) 


where fı and gı are the counterterms. 


C. Power spectrum 


Calculations can be simplified greatly in the spatially flat gauge (w = 0, E = 0), because 
the perturbed triad vanishes (0E? = 0) in this gauge. Here we fix the gauge after having 


9 
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put quantum corrections in Hamiltonian and having checked consistency. In contrast, in 
references " the authors fixed the gauge beforehand. We believe our treatment is more 


consistent. By solving the constraint equations, the perturbed lapse function and shift vector 


read 
| lag 1 
? 7 T ae 
1ā11 ; 
dm us LI (o dp — $6 fi)é c PPV o(p} -3H( --f)ó. (37) 


And the extrinsic curvature is 
ao Ki = —óVA(1 + f)o- 0,0! B, (38) 


where ak = H. 


In the spatially flat gauge, the total second order Hamiltonian becomes 


3pà 


-2 ca 
(2) gp epo. GODS un. " 
KH [a cl) 153 ;g gı 2h — 2f) 


2p? 


a ef 
+ EMV ao nenne Noe 


V 


ee E ab 
7 gı)ðp + S (1+ 95) Andparsy , (39) 


where Ó is also a correction function for the inverse volume and à? = 70. In this gauge, the 
dynamical inflaton perturbation dy coincides with the Sasaki-Mukhanov variable u = zt, 
with 


E - E + E — vo) SÉ f (40) 


Then, one can derive the Mukhanov equation 


n 


vi — (QA ell =: (41) 


where c, is the propagation speed of the perturbation. The solution of the above equation 


is |19 


H : 
u(k, T) ux om | + ikr - cp tA] ; 


Hew [r (2) sar (M)rrow). ^ vm 


10 
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where : 

F (2) = h - ct -[i«c D NW (43) 
where x = ovo(1 + 0/6)/3 + ao(5 — a /3)/2, C = BEC and ðo = 6(ko)/ao. The latter 
variable 6(ko) is constrained by the cosmic observational data bi). such as Cosmic Microwave 
Background and Large Scale Structures. For the specific inflationary models with a quadratic 
potential and o = 2, (ko) ~ O(10^?). In this paper we take ag ~ O(1077), i.e., the variable 
Au is of the order O(10~%). Moreover, Eq. (42) tells us that, on the one hand the inverse 
volume corrections become important for long wave modes with k < ko; on the other 
hand long wave modes cross horizon earlier than the short ones. It means that the inverse 
volume correction will leave more hints on large scales than small ones. These features are 
much different from those of the inflationary models with higher derivative terms such as 
K-inflation in Einstein gravity. 


Using the canonical quantization, we have 


CIE ai CT +07 CIE rag + C' (E, rah (44) 


Po 
where CIE. T) = u(k,7)/z(k). Then the two-point correlation functions of curvature pertur- 


bations can be calculated straight forwardly as 


(Cece) = UE ee (am) 96h + EL 1 + s) (45) 
kok) "" z2 FS T EENS ^J sOinv f » 
with 

Ys = vo(1 - 0/6) + a09/2e — x/(o + 1) (46) 


=> 


where the first and second terms in (46) comes from the z(k) factor in the gauge transfor- 
mations, and the third term attributes to the modifications of vacuum state (42). 


Finally, we can get the primordial power spectrum of curvature perturbations as 


H" au oh. (47) 


~ An?eg 
This result is in agreement with that in . When all the corrections vanish, this result 
is back to the case of the single field inflationary model in Einstein gravity 34). The 
primordial power spectrum of curvature perturbations and angular power spectrum are 
plotted in Fig. [I] where the dotted (purple), dashed (deep blue) and the solid (light blue) 
curves correspond to different values of the parameter C (C = 0,4 x 1074,3 x 10?) which 


was defined after Eq. (43). 


11 
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Pe(k)/P ko) 


(b) 


0.9 


5x10-^ 0.001 0.005 0.010 0.050 I 5 10 50 100 500 1000 
kMpc^! i 


FIG. 1: Primordial P;(k) (a) and angular C; (b) power spectra. Here, we set the slow-roll param- 
eters as the typical values eo = 79 = 0.01 and o = 2. In our calculations, the pivot wavenumber ko 


equals to 0.002Mpc?. 


IV. BISPECTRUM 


In this section, we firstly derive the third order Hamiltonian in the spatially flat gauge; 
then calculate the three-point functions of the primordial curvature perturbations; and fi- 


nally figure out the sizes and shapes of the bispectrum. 


A. "The third order Hamiltonian and In-In formalism 


We can get the corrected third order Hamiltonian 


FH = HOIN + HOIN) el, LN] AED, [BN] 


grav 
1 - 1 S 
=5 | d'eäianfl, +5 f d*zóNasQ., + f dx rap +052 +42) 
x x x 
n d Prd NH + 59 + O99) , (48) 


where the expressions for HP rw, 52), EE 8e ) are complicated, and we list them in Ap- 


pendix [A] 
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The perturbed third order diffeomorphism constraint is 


go] — 1 f 
Y 


a 1 ij E j l; c 
i drôN EE WEI + 40,6 Kj E! — 35 (BO. ELSE; 


— 0 KI5B° — Set — yôKjð ðE? + 9rd, (49) 
We calculate the non-Gaussianity in the interaction picture 
(Or) = (Us Cr) Uis (r, 70), Ui = ero e, (50) 
Up to the first order, we have 


(C(r) = -i f dr (I (7), 26s (7))- (91) 


TO 
B. Sizes and shapes 


Based on the second order anomaly free perturbative LQC theory, we combine Eq. (48) 


with Eq.(49) and arrive at the third order interaction Hamiltonian with some counterterms 


p 77 
MS = HO + go = | Pa| (1 — gig’ + E + g) aSpa 
p a / 3p 243 
+=(1 — g1)(0°B)dy'd.d~ — —(1 + 2f)H'd 


a 


D 
D d 1+ p? 8 
ECH +2fı— gı) ECH RÀ LE (oyhotó 
p = LE NE p Ae 
ura 2 = Qt fi- gue oy D Së Zeil (52) 


where the first three terms are in the leading order under the slow-roll approximation. In 
the following calculations, we only consider these terms. Comparing Eq. with those 
in 34], we could attribute two kinds of modifications in the interaction Hamiltonian to 
the inverse volume corrections. One comes from the modification of vacuum state (42), the 
other from the background dependent coefficients, such as (7,0,à,9g1,:--). However, the 
non-Gaussian signatures from the latter are contaminated greatly by the cosmic variance. 
Hence, we ignore the modifications in 7 etc., when we calculate the parts of three-point 
functions directly from the In-In formalism. In another words, we ignore the quantum 


anomaly behavior for the Hamiltonian in this work. In order to demonstrate the reasons 


13 
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more manifestly, we take the first term in as an example. 
] ^3 - g)5e" à 
= fas oP igi i: fox’ ae P Dër ( i 1) Jos oe" (53) 


where the first term in the second line appears in the usual form, while the second one 
contains An = Ann". Note diny depends on k, there is a corresponding k related to Aan, 
When we put these terms into the formalism and perform the time integral, we will 


obtain such corrections 
2 
ead Un kz, ks) +++: , K= ki + k+ ks, (54) 


where the dots denotes higher order corrections. From the above expression, we can see 
that such term peaks at k >> K. Since larger k gives smaller ópj, and quadruple is the 
lowest detectable mode in CMB, k corresponds to l = 2. So Eq.(64) peaks on the very 
low | < 2 region where the cosmic variance dominates over the signals (See Figure 1 and 
Figure 2 in LL Although the non-Gaussian features are presented on both large scale and 
small scale, above analysis implies that we can ignore the effects on small scales. Hence, 
we can ignore such terms in our following calculations. Once we ignore the corrections in 
the background dependent coefficients in Eq.(52), the form of the third order interaction 
Hamiltonian reduces to the usual one [30H34]. 

In order to eliminate the terms proportional to the linear equations, we need to do the 


field redefinition 


i E 
c 6-5 (1- fx) e et iso eno 
= Get OC + C207 (Geo Ge) . (55) 


Then, the interaction Hamiltonian can be reduced into the simple form 
HQ = [dopo es (56) 


According to ad , after a field redefinition of the schematic form ¢ = ¢, + AC? then 


the correlation function will contain two terms 


(C (23) 6 12) 6(23)) = (Ge) 6 (12) 65 03)) + 2A [(CG)6 (w2)) (C7) G23)) + cyclic] , (57) 
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where the first term can be computed by the In-In formalism and the second term comes 
from the field redefinition ¢ = C. + AC?. 
Let us firstly calculate the first term 


(Gr eG Ra) rl) = =i f eiis bl aol ben 


ror Ars fiho (eT) 
(58) 


where z(k) and C are defined in Eqs. (40) and (48), respectively. Here, we emphasize again 
that corrections from the background dependent coefficients are neglected in the above 
expressions. However, in the next calculations of the parts from field redefinition, such 
corrections should be taken into account for consistence. Because the background dependent 
coefficients (i /?1|, ,---) take the value at the moment when the corresponding mode crosses 
horizon (T, ~ k~'), they contain corrections such as (ko/k;)?. Of course, these terms also 
peak at the points where kg > k;, they become important, particularly, in the squeezed 
triangle limit (kj < ko, ks). 


The contributions from field redefinition can be decomposed into two parts, one is 


E (59) 


Hi 
(27)?6( k;) MEME 
8 E SE EIS 


and the other is 

(2)°5(}— &) en Zi X aos] (60) 
where H — 8 and the overdot stands for derivative with respect to the cosmic time. The 
C,(k;), Co(k;) terms can be read from Eq. (55) by substituting 7 with k; '. 

In summary, we conclude that the forms of interaction Hamiltonian are exactly the same 
as the usual one 34], however, the inverse volume corrections ðiny will make some con- 
tribution to the bispectrum. There are mainly two sources, one is the modifications to the 
standard Bunch-Davies vacuum, the other comes from the z(k) factor in the gauge trans- 
formation Q(k) — u(k)/z(k). Furthermore, we can expand the above results in terms of 


Oinv = Óg(ko/k;)" and obtain 
(2x)58(V ki) E angle (Ki, Ka, k3) + Fi (kı, ko, ka) + Zb, ka, ball. (61) 
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where 


Zar Hl P? p 
Jum je = _ re d (3€o — 270) ke + en k; ki j Bé Doing kik en" 
i 3 K 
4 ITA QD) - 


is the usual leading term and 


r= D ol us (63) 


(2m)4P? afko y a, (bah 
F = [TE o] i | 204(¢o =) + wheo| 2 k; (=) ZE » SC (2 | 
(64) 


are the inverse volume correction terms, and here the relevant coefficients are 


wy, = 20—C,, 65 


Wy = Ye + C,-2C, 66 


(65) 
(66) 
US = yq C; —2C, (67) 
(68) 


c 7 OQ 
z= — compa a 
C Au |^ L + ) + Gi 68 
PESE (69) 
2(c +1) 


Particularly, we figure out that the corrections from Bum z()] terms in the denomi- 
nators in Eqs.(68), (59) and (60) are absorbed into F, shape, while all other corrections are 
collected in Fy. In the next section, we will find F; shape provides an unique signal from 
LQC mechanism, and more importantly, this signal is independent of the inflationary mod- 
els, because MS 214 terms always appear in the gauge transformations (40). Namely 
it is an universal signal in the LQC scenario. 

The single shape is the usual one, while JF; and F> arise only in LQC scenarios. Sizes 
of the two new parts are proportional to parameters (w; ,:::), which are in the order of 
O(107?, 1072). That is to say that, these new non-Gaussian features from LQC are smaller 
than those of usual inflationary models in Einstein gravity by a factor at least 10 ?. Although 
this does be a tiny number, considering these features sourced by quantum effect, this factor 


is not small as initially expected. Furthermore, as stated above, because the inverse volume 
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mi, 


0.5 


1.0 X3 1.0 X3 
FIG. 2: Single (a) and Fı (b) shapes. The z-axis is 2533 (x2, 13) / P2, here we set z; — 1. The 


slow-roll parameters take the typical values eg = ro = 0.01. Here and hereafter we set o = 2. 


corrections in the interaction Hamiltonian can be neglected, i.e. we could use the usual 
Hamiltonian directly in the In-In formalism by substituting the Bunch-Davies vacuum state 
with the one given in Eq.(42), we argue that the expectations on the sizes of bispectrum 
should hold for any inflationary models in LQC scenario. 

In Fig. 2] (a) (b) and Fig. Bl(c), we plot z1z5x3J (zi, x2, £3)/Pê, with z; = ki/k. 
The difference between shapes Fı and F is plotted in Fig. (d). We can see that all 
the three shapes peak at the squeezed limit (xə = 1, r3 = 0), however, the substructures 
are different among them. Compared to the single shape, Fı shape upraises at the corner 
(z2 = 0,23 = 0.5), while F> flattens at the same point. From Fig. [3l (d), we can also see 


that Fə peaks more dramatically in the squeezed corner than F4. 


V. SHAPE CORRELATIONS 


In the Figures [2] and B| all the three shapes looks similar. In order to figure out the 
differences among shapes more quantitatively, we need to calculate the correlations between 


them. 
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(d) E. 


1.0 X3 


FIG. 3: Fə (c) shape and Fy — Fı (d). The coefficients (wi ,w5 ,w$) in Fi, Fo are set unity. 
Firstly, we define a 3D shape function 40 
1 
S(ki, kə, ka) = Uto) Fla, kə, k3) ; (70) 


where N is a normalization factor which will not affect the following calculations. 


Then, we construct the products of two shape functions 
asi — | Ss. ko, ka) S' (kn, ka, jon a a) (71) 
Vk 


where w(ky, ko, k3) is a weight function and V; is the integration domain constrained by the 
triangle inequality. 

Finally, we arrive at the 3D shape function correlator 
F(S, S") 


/F(S, S)F(S', S") - 


This quantity describe the cross correlations between two different shapes. 


C(S, S") = (72) 


Furthermore, for a large series of well-motivated shapes, the above descriptions can 
be simplified. One can define the distance described by k from the origin of (kı, ko, k3)- 


momentum space to the particular triangle slice which is perpendicular to (1, 1, 1) direction, 


k E (ky + ko + k3) ; (73) 


N| = 
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then, we introduce another two new variables 


iom sk tats), (75) 
cg = SEH (76) 


In the domain constrained by the triangle inequality, 0 € k € oo, 0 € 8 € 1 and -(1— 8) € 
o < 1— B. For the classes of models with homogeneous shape, which means that the powers 
of wavenumber in shapes are homogeneous, the k dependence in the 3D shape function can 


be separated 
Sky, ka, k3) = f(k) (o, EI. dV, = dkidkadks = k?dkdadf . (77) 


In fact, for the models considered in and here, f(k) « const. Hence, we can focus 
on the 2D shape function Z (a, 8) and the integral over k cancels when one calculates the 
shape function correlators. 


Further, one can introduce 
a=1-2, B=ay/3, O<2,y<1), (78) 
to square the integration regime. By using these variables, the integral measurement becomes 
dad’ = rdrdy . (79) 


From the above expression, we can read the weight function w(x, ul = x. This choice works 
well for all the shapes mentioned in |40H42], however, for our new shapes Fı and F> the 
correlation matrices Cm, defined in (84), suffers from divergence. In our calculation we 


therefore use a new variable £? = x 
dad = @dédy , (80) 


to eliminate such divergence. Using variables (£, y), we can decompose the shape Z (£, y) on 
any triangle slice with analogous radial polynomials R,,,(€) and shifted Legendre polynomials 
P.) 

S (E) — M en Bal El P (y) ; (81) 
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where the first few R,,(€) eigenfunctions are 


Ro=V2, Ry, = V4(-24+36), Ro = V6(3— 126 + 10€?) , 
Hs X 8(-4d 4 B06 — 6022 +358) ee (82) 


And P,(y) eigenfunctions are 


Po=l, P=V3(-1+2y), Ps V5(1-6y+ 6y’) , 
P, = (1 + 12y — 30y? + 204°) ,---. (83) 


Thus, one can define the correlation matrix Cmn 
1 1 7 
com c f AE f du (E) RA (OP. QD. (84) 
0 0 


The Cmn matrices for local, equilateral, single field model, JF; and F, are listed in (85) and 
Fig. Æ| The explicit definitions of 2D shape function .7/(£, y) can be found in Appendix [B] 


1.00 —0.14 0.03 0.00 1.00 0.40 —0.12 0.01 1.00 —0.07 0.01 0.00 
0.38 —0.07 0.02 0.00 0.68 0.27 —0.09 0.01 0.43 —0.02 0.01 0.00 
0.04 —0.01 0.01 0.00 | | 0.21 0.07 —0.03 0.00 | | 0.07 0.00 0.00 0.00 | - 
0.02 0.00 0.00 0.00 0.01 0.00 —0.01 0.00 0.02 0.00 0.00 0.00 
1.00 —0.19 0.04 0.00 1.00 —0.12 0.24 0.00 
-0.26 009 —0.01 0.00 | 0.30 —0.03 0.01 0.00 fas) 
0.44 —0.10 0.02 0.00 0.07 —0.01 0.01 0.00 
—0.32 0.08 —0.01 0.00 0.02 0.00 0.00 0.00 


From Fig. [4] we can see that JF, term differs from others explicitly, while F> is almost 
indistinguishable with the local form visually. 


Armed with the above results, one can calculate 2D shape correlator 


oz / 
CS, PN = — IIZI , (86) 
F LS un un 
where the product is defined through 
Et, I=] FEN E dY ënn (87) 


Sk 


m,n 
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Ry) 


ce oe ee 


1 | 0.826] 0.996 | 0.722 | 0.996 
1 | 0.874 | 0.535 | 0.826 
1  |0.707 | 0.992 


TABLE I: 2D shape correlators. 


The numerical results of the 2D correlators are listed in Tablet |l] from which we can find 
that the correlations of shape F, are low with all other four shapes, while shape Za possesses 
high correlations with local (particularly), equilateral as well as single form. That is to say 
that Fə is almost indistinguishable with local form, while JF; does be the unique signal of 
LQC. As already argued in the previous section, Ei is an universal shape in LQC scenario, 
so we can identify this shape as a new window for LQC scenario. 


Finally, let us estimate the parameter fnt. Here we focus on the two new shapes F12, from 
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Tab. [I]we can see that Fo is highly correlated with the local form, while F, is less correlated 


with them. The contributions to fX?! from these two shapes can be easily estimated as 


AJNE ~ 9 X G (Foca Fria) X fx" (single) ~ O(10 7) fr^ (single) , (88) 


where w represents in short parameters (w1, w5, w5) in and (64), whose typical values 
are of order O(107?). Because f} (single) in the usual case is of order ~ eg |43], the 
contributions from inverse volume corrections is completely negligible. However, as argued 
before, our results should be robust for other inflationary models in LQC scenarios, especially 
those with large non-Gaussianities, such as K-inflation, DBI inflation etc. So, one can 
anticipate that in those models with large non-Gaussianities the features from inverse volume 


corrections in LQC scenarios might be observed. 


VI. CONCLUSIONS 


In this paper we investigated the contributions to the cosmic primordial scalar bispec- 
trum from the inverse volume corrections in LQC scenarios. We derived the interaction 
Hamiltonian, however, we found that the new interactions contribute greatly to the modes 
with kı + kə + ka « k. Because the scales corresponding to k is very large, here we take 
k ~ 0.00014Mpc ! corresponding to quadruple mode l = 2, it means that the three wave 
lengthes in the bispectrum are all on the super-horizon scales. On so large scales the cosmic 
variance usually dominates over the signals, we hence neglected these new interactions in 
our calculations. That is to say, the interaction Hamiltonian we used is the same form as 
the usual one for the single field inflation models in Einstein gravity. This greatly simplifies 
our calculations, and more importantly, makes our results robust, i.e., our results should 
hold for other inflationary models in LQC scenarios. 

Although the Hamiltonian shares the same forms as the one in Einstein gravity, the inverse 
volume corrections diny X (ko/ k)^ still contribute to the bispectrum. Roughly speaking, there 
are two aspects, one comes from the deviations from the standard Bunch-Davies vacuum; 
the other attributes to the non-trivial gauge transformations ¢(k) = u(k)/z(k) from the 
spatially flat gauge to the observable curvature perturbations. Consequently, we obtained 
the three-point functions of the gauge invariant curvature perturbations. We found that, 


except for the usual single component in slow-roll inflation models, two new shapes arise due 
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to the corrections, namely 75. Furthermore, we performed a careful analysis on the new 
shapes. We found that, the whole profiles for all the three shapes (single, Fı,2) are visually 
similar, i.e. they peak at the squeezed limit. However, the substructures among them are 
different. Compared to the single shape, Fı shape upraises at another corner (See. Fig. [2] 
and[3), while F, flattens at the same point, and F, peaks more dramatically in the squeezed 
corner than Ji. In addition, we investigated the correlations among five shapes, including 
local, equilateral, single and Pia, The results show that Fə is highly correlated with the 
local type, while F; is less. It means that the latter can provide a new window for probing 
the loop quantum mechanisms using cosmic primordial bispectrum information. Finally, we 
estimated the order of observable parameter A flee! ~ O(1073) x fyec!(inflation) from the 
inverse volume corrections. 

The non-Gaussianity from the inverse volume corrections in LQC scenarios is tiny and still 
undetectable currently, however, considering they are generated by the quantum effect, which 
is naively expected of the order O(GUT/Planck)? ~ O(10 ?)^, our finding becomes non- 
trivial. Especially, the results obtained in this work should be generalized directly to other 
inflationary models with large non-Gaussianities, in which the inverse volume corrections also 
becomes large therein. In addition, in this paper we only investigated the non-Gaussianities 
from the inverse volume corrections, while ignored those from the holonomy corrections, 
in which the sound speed of scalar perturbations are typically changed. Besides, we only 
investigated in this work the bispectrum from scalar modes, left those from tensor mode 


unexplored. These topics are worth investigating further. 
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Appendix A: Interacting Hamiltonian 


In this appendix, we derive the perturbed Hamiltonian density and the perturbed dif- 
feomorphism constraint up to the third order. According to , the classical Hamiltonian 


includes two parts 
JEN] = Hren [N] F GER [N] (A1) 


The gravitational Hamiltonian can be expressed in terms of the extrinsic curvature 


ASIN] = di deN% = d para EIER 20 rage (rere — K"K7)|, (A2) 
grav 2 2 6 | det E] Cd mn\~e~d c d/|? 
and the matter part of the Hamiltonian is 
desch) EC f PEN (Hr F Ny F Ho); (A3) 
S 
where 
ar? _ Eè Gecke s 
Ör = >, Ny = — = y | det E|V (v A4 
24/ | det E| 2,/| det E| ' | Va 


The diffeomorphism constraint is 
d eas [IN] i [! deN’ (a Al — 8, AJ) Eh — Ai Ia, (A5) 
and the matter contribution is 
f ashes iN | = [ esae. (A6) 


In order to obtain the perturbed Hamiltonian density and the perturbed diffeomorphism 
constraint up to the third order, we need the following two relations. We expand (det E)? 


and (det E)3 to the third order as follows. 


1 3 ie LL -— 
E)?=p?| | 14+ 68Et-G LED? — —HÓóiójoERÓES + — (8*8 Eg)? 
(det ) p | + 2p a i + ap a W App? b~a i j i3 Bp c k) 
1 c Zei a 1 i j a c 
- sr ; (056 Ez) (8628 E26 E) + — EAE 
+ ES Sils ES ES6ES + (A7) 
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(det E) + =p] 1- iE; + gaS) + fg - —l _(si5 23 


48p? 


(656 ES) (65020 E 6E) — EN Op LO ESO EO Ey 


1 
8p3 12p3 ° 


= — aiio S Et GEIGER + 
127? 
Thus the third order gravitational Hamiltonian density can be written as 


SOT SP EAS EOD 


where 


ELT Lët, 2e] 8 ES Ej [(det E) 2] 8,9T 4 


i diu [(det EI" 2] 8,6T4 + 22 ESE (det E) 3] 0,61 


= onser SES(O.0-6 Et) + 2 (0.026 Eg) 
E 
— 1 se(óm8E? J(0,8,8E7) + EE 
Ap? 2p3 


c md 
59 = Gi Ej Ek el rr (3) 


2 det E Emntctd 
mg , ESES ` 
EET Re 
e e 
rel Et Ed[(det E" 5] i, 68" 


STT” 


" 1 SS 
m 10 ERO (8,0 E) (0,0 Ez) - 5519 9 E(OAEL (SED), 
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= — SO ERO E OK) + 


EE, pem. ft PREIS SKPK! 
= —€ Emn c €; Emn c 
ydet E S v det E S 
Be 4 ESSE? ` 
jk J k i m n jk ~j k d m n 
€ meer S s DA e ap te òK; 


—e)"S E56 Ej [(det rely 


-ef SES Ei (det E) SI ei, OK K} — SEER (det E 


i Mi 


LS Ah" — d  ES5EQ(det E 


) Zu 
—e* BSS Ef] (det E)-3] 
) 2 


—e* 45 E(det EI SI es — ef" BSS Ef (det E)-3 


—e)" E Ej [(det E) 3]Be SK” R? — c^ E Bal (det E) joe. 


=P E Ed((det E" ll KOK? 
2k 2k 
us — JO ESO EROKE — 


ns 


No Bio 


EE óECOKT OK) + 


2 1006 0 ESO ERO 


E A 
F LAAEIT — 


(éi KLYPE GE? — l gegi OK OKO SE? 


l"n"m 


p2 


j vd c n sl a 1 k? 
a EE 
p2 


"g Et)’ (618K?) — 


2k 
+ = 06E 
p? 


a 
2 


(60E) (K7) — 


SÉ S SELE?) SES) — (6 


EE 


2p? 
(ESEE + 39 (alan) tate Ets BS) 
sp: 4p? 
k Hi 
"O10 535 B86 BX SES + — di 5158 SELS ELS EC 
NES 2p2 


The matter hamiltonian can be expressed as 


$9 


1 1 1 
= 561° [(det E) SI + zón[(det E) 3]? + 
(Ab 
8p? 4p? 2p 
(GOLS EOE) (HSE!) SOR SIS ESO EO EE 


-Sins Ee + Tôn] 
Ap2 p? 
l (Ab 
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,ET K? — d" Ez Ed[(det E)- 3] 


dà OK™5K™ 
rue On 
cue ee? 


Senis Hg 
KrK? 


2 f 
— OK POESK) 
p2 


Er 


SKILE? 


(025 E; 


1910 EF 6E?) (545K?) 


(A12) 


7? [(det E)-3]9) 
sisi RES z? 


3 
o IF B25 E* E; 


Do | (A13) 
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09 EPO E20,5p0,59 ELSE! ð pð ` OV EE? 


59 si) det E)-3 (09. 60,6 
y 2v det E 2y det E 2 K r*j get? 
ai b k c Sab 
_ 6M OE; EE | O00 E06 VIVUM (A14) 
p? Ap? 
1 1 1 2 1 8 
HD = gj Vdet EV, (p)ov?  z[(det E)5]? V. (z)ó" + [(det E)?]°V (9) 
178) 
+ [am| vue 
13 8 E tee 7 ô SEL)? 
= SDV oolp)" ga qP OE Vaelëlä e [oem 
8p2 
6,020 E06 E" E a | (ESEC) — (OF Ef) (5,025 E76 E?) 
m Api Vo(9)óo "m p 4803 a 8p? 
OOF OHSOR OES — 01036 O E26 EO EG 7 
De ps a (a15) 


Combined with the results obtained by ld. we get the perturbed Hamiltonian as follows. 


AMAN] = PASIÓN] + DER N], (A16) 
where 
MRIN = 5 | Pav He, (A17) 
and 
HREN] = 5 | PENIG (A18) 
and the matter Hamiltonian reads 
Hall = | PeR [59 +99 +20), (A19) 
Hiel N] = | d'ssN [n +59 +90]. (A20) 


Here we have ignored the high order correction terms caused by the inverse volume, such 
as a?) HO, because in the in-in formulism, these terms do not contribute to the non- 
Gaussianity. 


The diffeomorphism constraint up to the third order is 


1 a 7 1. ; 
2 [SN] :— = H d?zóN^ EE zer, CAERE; 
: Lo c ; 
— O,5KL6E? — zer, DEROSE; — y6K1à E), (A21) 
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29 aló N?) = d d*xóN* 60,6. (A22) 


Appendix B: Definitions of shape functions 
For all the shapes considered in this paper f(k) — const., so we have 
S(Ki, ko, k3) = Z (€, y) . (B1) 


For the local shape 
ke + k3 + k3 


ocal — B2 
Stocal Bk koks (B2) 
for the equilateral shape 
(ky + ko — k3)(ko + ka — kı) (k3 + kı — kp) 
equil — —= -—-————— —' B3 
im SS ge 
for the single shape 
EIS? — 210) 25 k? + €o b kik? + Ben Xi EK 
Ssingle = —— U (B4) 


kı ko ks 
ko N^ ko M" ko \° 
S FEE TUE ue Dan e» B5 
an (2) dës +(R) | "SS E? 


and for the Fj shape 
Ai — no) + wheo| 3 k? (Ko/ ki)” + waco 3 2,4; kzka (ko/ ki) 
Srn = — ————————————— . (B6) 
ky Ko ks 


for the Fı shape 
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